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1 Introduction 

It is a well known result that when electrons move inside semiconductors, they experiment 
the potential barriers due to the heterostructures of the material. With basics quantum 
mechanics, by solving the Schroedinger equation we can compute transmission and reflec- 
tion coefficients, making possible the understanding of the resonance structures of each case 
of potentials studied. This is one of the aspects of theoretical condensed matter extensively 
studied in the 90's [1-3]. Due to the improvement of numerical computational methods, 
a few more difficult problems without exact solutions were attacked. The source of those 
difficulties is the complexity of the several kinds of potentials that should be treated in 
the Schroedinger equation. In particular, a numerical method was developed in order to 
solve the problem of potential barriers for any case [4]. In order to obtain these results 
numerically, they made use of a common tool in condensed matter systems: The trans- 
fer matrix. In a few words, the program computes transmission coefficients for the given 
potential profile. 

The fact is that, in the last years, the models in physics of extra dimensions brought 
back to the arena the analysis of the Schroedinger equation with a potential generated 
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by thin membranes. In these models, the interaction between the several fields and the 
membrane is an important aspect to be understood. For example, localization of fields and 
the structure of resonances can tell how the the universe inside the membrane behaves. 
In that sense, particles moving through the extra dimensions resembles the wave packets 
of electrons hitting the potential barriers inside semiconductors. The core idea of extra 
dimensional models is to consider the four-dimensional universe as a hypersurface embedded 
in a multidimensional manifold. In particular, the standard model presents interesting 
topics to study such as the hierarchy problem, and the cosmological constant problem that 
can be treated by the above-mentioned scenarios. For example, the Randall-Sundrum (RS) 
model [5, 6] provides a possible solution to the hierarchy problem and show how gravity is 
trapped to a membrane. 

For compact extra dimensions, as in Kaluza-Klein models, a massless field in five 
dimensions gives origin to a massive one in four dimensions. The spectrum of masses is 
discrete and they are determined by the radius of the extra dimension. In RS models the 
extra dimension is infinity. In these models, generally the mass spectrum is determined 
by a Schroedinger like equation and the spectrum is not guaranteed to be discrete. In 
these models it is important to verify if electromagnetism and gravity are reproduced in 
first approximation inside the brane. This is true if the zero mode, or massless mode, 
is localized in the membrane. In these scenarios, some facts about localization of fields 
are known: The scalar field (0— form) and the gravity field are localizable, but the vector 
gauge fields (1— form) is not. The reason why this happens is that, in four dimensions, the 
vector field theory is conformal and all information coming from warp factors drops out, 
necessarily rendering a non-normalizable four dimensional effective action. Beyond this the 
left hand of the fermion field is also localizable. 

It is important to note that the presence of one more extra dimension (D = 5) provides 
the existence of many antisymmetric fields, namely the two, three, four and five forms. 
However, the only relevant ones for our membrane are the two and the three form. This 
is due to the fact that when the number of dimensions increase, also increases the gauge 
freedom. This can be used to cancel the dynamics of the field in the visible brane. From 
the physical viewpoint, they are of great interest because they may have the status of 
fields describing particles other than the usual ones. As an example we can cite the space- 
time torsion [7] and the axion field [8, 9] that have separated descriptions by the two-form. 
Besides this. String Theory shows the naturalness of higher rank tensor fields in its spectrum 
[10, 11]. Other applications of these kind of fields have been made showing its relation with 
the AdS/CFT conjecture [12]. 

For these fields it is also a known fact that their zero mode are not localized. The mass 
spectrum of the two and three form have been studied, for example, in Refs. [7] and [13]. 
Posteriorly, the coupling between the two and three forms with the dilaton was studied, in 
different contexts, in [14-16]. The study of this kind of coupling, inspired in string theory, 
is important in order to produce a process that, in principle, could be seen in LHC. This 
is a Drell-Yang process in which a pair quark-antiquark can give rise to a three (two)-form 
field, mediated by a dilaton. 

In this approach, we make use of soliton-like solutions: these are studied with increasing 
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interest in physics, not only in Condensed Matter, as in Particle Physics and Cosmology. 
In brane models, they are used as mechanism of field localization, avoiding the appearance 
of the troublesome infinities. Several kinds of defects in brane scenarios are considered 
in the literature [17-20]. As an example, in a recent paper, a model is considered for 
coupling fermions to brane and/or antibrane modeled by a kink antikink system [21]. The 
localization of fields in a framework that consider the brane as a kink have been studied 
for example in [22-26]. 

A proposal to solve the gauge field localization problem was given by Kehagias and 
Tamvakis [27]. They have shown that the coupling between the dilaton and the vector gauge 
field produces the desired localization. In this work they also consider the localization of 
gravity in the background with and without the dilaton. In analogy with their work, the 
present authors have also considered the issue of localization and resonances of a three- 
form field in [28] and in two separate papers, g— form fields in a scenario of a p— brane with 
co-dimension one and two [29, 30]. As has been pointed out in [29], the coupling here is 
needed for a different reason. In that work, it was shown that the some form fields without 
a dilaton do not have its massless modes localized. When we introduce a dilaton, this mode 
can be localized. 

Therefore we are left with two phenomenological setups: with and without the dilaton 
field. In the case of form fields, only the zero form do not need of the dilaton coupling. 
When we consider the other fields, namely the gravity and fermion fields, the dilaton is 
not needed for localization. In the recent literature, generally the dilaton is considered in 
the form field case and not in the gravity and fermion cases. Therefore a more complete 
analyses should consider all the fields in both setups. 

In all the models cited above a common feature is that the massive modes are not 
normalizable and therefore we do not have a well defined effective action. Despite of 
this there is the important possibility of appearance of resonances. As said before the 
mass spectrum is determined by a Schroedinger like equation. The potential term in this 
equation will generally have a volcano shape and therefore we can ask about the possibility 
of massive modes, besides living in the extra dimension, have a peak of probability to 
be found at the location of the membrane. This analysis have been done extensively in 
the literature [31-46] . For example, in [45], the authors consider a the vector field and 
found that, after normalization of the the wave function in a truncated region, there is a 
resonance very close to m = 0, and this indicates that the lighter is the mass, the bigger is 
the probability of interaction with the membrane. 

From our viewpoint the wave function has an oscillatory part and cannot be normal- 
ized. In order to study this, the authors define a relative probability, but for us it is not 
clear how it solves the problem. Another very important point is the truncation of the in- 
tegration region. This is equivalent to have a potential that do not falls to zero at infinity, 
therefore it is natural to have a peak of resonance very close to m = 0. This could lead to 
the wrong interpretation that very light modes can interact with the membrane. From our 
analysis, it is not true that the lighter are the massive modes, the bigger is the probability 
of finding them inside the membrane. In fact, using our method it has been show [28, 29] 
that heavy modes also can interact with the membrane differently from what is discussed 
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in [46]. In fact we have found a very rich structure of resonances. 

In order to analyse resonances, we must compute transmission coefficients (T). This 
gives a clearer and cleaner physical interpretation about what happens to a free wave that 
interact with the membrane. The idea of the existence of a resonant mode is that for a given 
mass the transmitted and reflected oscillatory modes are in phase inside the membrane, i. 
e., the transmission coefficient has a peak at this mass value. That means the amplitude 
of the wave-function has a maximum value at z = and the probability to find this KK 
mode inside the membrane is higher. 

Therefore the main subject of this article is to give in detail the numerical method 
used recently to study aspects of the form fields in a scenario of extra dimensions. As 
said before, it is lacking in literature the study of resonances for form fields without the 
dilaton and for gravity and fermion fields with the dilaton. Therefore as an application of 
the above method we compute the transmission coefficients of Gravity, Fermion and form 
fields in both setups: with and without the dilaton coupling. The paper is organized as 
follows. Section two is devoted to discuss a solution of Einstein's equation with a source 
given by a kink with and without the presence of the dilaton. In section two we give the 
general prescription to reach the Schroedinger equation used to analyses the problem. We 
also present in detail the numerical steps to compute the transmission coefficients with the 
transfer matrix. In Section four we analyses resonances of the gravity field. In section 
five, six and seven, we present the resonance structure for 0,1 and 2 forms respectively. In 
section eight we present the same study for the fermion fields. At the end, we discuss the 
conclusions and perspectives. 

2 The Kink as a Membrane 

We start our analysis by studying the space-time background. It is well known that vector 
gauge fields in these kind of scenarios are not localizable: in four dimensions the gauge 
vector field theory is conformal and all information coming from warp factors drops out 
necessarily rendering a non-normalizable four dimensional effective action. However, in 
the work of Kehagias and Tamvakis [27], it is shown that the coupling between the dilaton 
and the vector gauge field produces localization of the later. On the other hand scalar and 
fermion fields do not need this coupling in order to produce localized zero modes. Therefore 
we must consider in this section both backgrounds: with and without the dilaton coupling. 

2.1 The Background with the Dilaton Coupling 

In this section we obtain a solution of the equations of motion for the gravitational field 
in the background of the dilaton and the membrane. For such, we introduce the following 
action, similar of the one used in [27]: 



where G is the metric determinant and R is the Ricci scalar. Note that we are working with 
a model containing two real scalar fields. The field cj) plays the role of membrane generator 
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of the model while the field vr represents the dilaton. The potential function depends on 
both scalar fields. It is assumed the following ansatz for the space-time metric: 



ds2 = e'^^'^y^7]^^dxf'dx'' + e'^^^y^dy^, (2.2) 

where rjf^i, = diag{—l, 1,1,1) is the metric of the p-brane, y is the co-dimension coordinate. 
As usual, capital Latin indexes represent the coordinates in the bulk and Greek indexes, 
those on the brane. The equations of motion are given by 

1(<^')2 + l(7r')2 _ e2Bfo)y(^^^) ^ 2AM'\A'f, (2.3) 

^{(P'f + f + e^^^y^V{(t), vr) = -12M3^" - 2AM^{A f + 12M^J^B', (2.4) 

+ (4yl' - B')<p' = d^V, (2.5) 

and 

vr" + {AA' - B'y = d^V. (2.6) 

In order to solve this system, we use the so-called super-potential function W{(j)), defined 
by 0' = following the approach of Kehagias and Tamvakis [27]. The particular solution 
regarded follows from choosing the potential V{(j), vr) and super-potential W{(j)) as 

and ^ 

W{ck) = vact> {^-l^. (2.8) 

where a and v are parameters that adjust the dimensionality. As pointed out in [27], 
this potential give us the desired soliton-like solution. In this way it is easy to obtain first 
order differential equations whose solutions are solutions of the equations of motion above, 
namely 

TT = -\f^APA, (2.9) 

At: , , 

B = ^ = ; , 2.10 

, W 

A' = (2.11) 

12M3 ^ ^ 

The solutions for these new set of equations are the following: 

(t>{y) = vtaiih{ay), (2.12) 

Aiy) = (41ncosh(ay) + tanh2(ay)) (2.13) 

All of these results are going to be used for analyses resonances of several fields with 
v^/72M'^ = 1. 
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2.2 The Background without the Dilaton CoupUng 

Here we look for the background solution without the dilaton coupling. Wc could repeat 
all the steps of the last section to arrive at the final solution. We can also take a short way 
and from our previous solution we obtain the desired result. By just setting i? = vr = we 
can get the final answer. Therefore we get for the metric 

ds^ = e^^^y^rj^^dxf'dx'' + dy\ (2.14) 

and for the resulting potential 

This solution agrees with the one obtained previously. We must note that the solution for 
A is left unchanged and therefore the contribution of the dilaton is to change the equations 
of motion. Therefore wc end here the analysis of this background and in the next section we 
present the general set to be used throughout this paper in order to look for the spectrum 
of resonances of fields. 

3 The Numerical Method 

In this section we give the general prescription to be followed in all the cases to be studied 
in this work. First we must show how to reach the equation that drives the massive modes 
and argue that it is a Schroedinger-like equation. Next we show in detail the transfer 
matrix method to compute transmission coefficients of an arbitrary potential profile. 

3.1 The Schroedinger-like Equation 

The first step is to find the equation of motion for the field in five dimensions, which 
generally takes the form 

6$(x,y) = 0, (3.1) 

where O is a differential operator em five dimensions. The next step is to separate the 
differential operator in the brane (O^d) and extra dimension (Oy): 

O = 04d + Oy. (3.2) 

After this we perform a separation of variables in the field ^{x,y) = il){y)(f>{x). It is 
important to note that when doing this all the information about the polarizations will be 
in (f) and ^ will be a scalar function. After these considerations we get 

Oid't'ix) _ Oytpiy) 
(f){x) tp{y) 

and as the Ihs depends only on x and the rhs on y they both must be equal to a constant 
that we call m^. The motivations for this constant is obvious after we write the equations 
of motions as 

OAd(t>{x) = -W?(l){x) 
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and 

The solutions of the second equation above therefore give us the aUowed masses in the 
visible brane. This equation in all cases will have the form 

+ P\y)jy + V{y)^ i;iy) = m^Q{y)i;iy). (3.3) 

The best way to study this equation is to transform it in a Schroedinger type equation 

+ U{z)]ij{z) = m^i;{z), (3.4) 



through the transformations 



dz — 

:r = f(y)^ 4^{y) = mm, (3.5) 

ay 



with 

fiy) = VQ(JJ), ^(y) = exp(P(y)/2)Q(y)-l/^ (3.6) 

and 

U{z) = V{y)/f + {P'{y)^'{y) - n"{y)) /^f (3.7) 

where the prime is a derivative with respect to y. In order to obtain a meaningful effective 
action we also need the integration in the extra dimension be finite. After the above 
transformation of coordinates, the interesting thing is that this finiteness is translated in 
the condition 

ip dz = finite 

which is like the square integrable condition of quantum mechanics. Therefore our problem 
can be addressed as if we were dealing with an one dimension quantum mechanical problem 
with potential given by (3.7). The study of the localization of massless modes (m^ = 0) is 
very simple and it has been widely studied in the literature. An obvious solution to the 
massless case is the constant mode ip = ipQ. With this the condition for localization reduces 
to 

/ Qdy = finite. 

The study of the massive modes must be done more carefully. This happens because when 
we go to specific cases we generally get that the potential has the asymptotic behaviour 

lim U{z) = 0. 

z— >iboo 

Therefore we do not have a discrete spectrum for m > 0. Beyond this, it is a known 
fact from quantum mechanics that any solution for positive m must posses a oscillatory 
contribution, and therefore the wave function is not normalizable. The conclusion is that 
the only possible localized mode is the massless one. This means that this is the unique 
mode living inside the membrane. 
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3.2 The Transfer Matrix 



Here we must give the details of the program used to compute the transmission coefficients 
by transfer matrix. First of ah, we consider as an example the simple case given by the 
double well barrier as shown in the Fig.l. 



U(z) 



Figure 1. The double potential barrier. 

The solution of the Schroedinger equation for each region is given by 

Mz) = ^ie""'' + Bie-'''^', h = ^/2{E-U^), i = 1,2,3,4,5. 

and is well known from quantum mechanics that it exhibit the following resonance profile 
as shown in Fig. 2: 

log(T) 




Figure 2. Resonance in the double potential barrier. 

The idea therefore is to apply the results above to a more general case. We are 
interested in apply that idea for an equation like 



(3.8) 



with potential given by the Fig. 3. 

In order to solve this, some authors have tried to find numerically relative probabilities 
defined as follows: 

As said before, the opinion of the present authors is that the above prescription hides 
the physics of the problem and the right approach should use plane waves, looking for 
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Figure 3. General potential with parity symmetry with lim2_j.-|-oo U{z) = 0. 

transmission coefficients. For this purpose the potential can be approximated by a series 
of barriers showed in the Fig. 4: 




Figure 4. Potential approximation by barrier. 

The Schroedinger equation must be solved for each interval Zi-i < z < Zi, where we 
have approximate the potential by 



U{z) = U {zi-i) = Ui-i, Zi-i = {zi + Zi_i)/2. 



with graphic given by the Fig. 5 
The solution in this interval is 



and the continuity of the tpi^i{z) and Tp'^_-^{z) at z = Zi give us 



Ai 



Mi 



A. 



Bi- 



(3.10) 



(3.11) 



(3.12) 
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Zi Zi-f-1 



Figure 5. The multistep regions. 



In the above equation we have that 

and performing this procedure iterativcly wc reach 

where, 

M = MnMn-i---M2Mi. (3.15) 
The transmission coefficient is therefore given by 

T = l/|M22p (3.16) 

and this expression is what we must to compute as a function of A, which in our case 
is m^. In order to obtain the numerical resonance values we choose the Zmax to satisfy 
U{zmax) ~ lO^'^ and m? runs from Umin = U{zmax) to Umax (the maximum potential 
value). We divide 2zmax 

by 10"^ or 10^ such that we have 10^ + 1 or 10^ + 1 transfer 

matrices. 



4 The Gravity Resonances 

As said in the introduction we are not going to study the localization of zero modes of 
fields. This have been done extensively in the literature. Here we must concentrate in the 
possibility of the existence of resonant modes. We consider the two backgrounds found in 
the last section. 
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In order to compute resonant massive modes we just need the Schroedinger like equa- 
tion that is obtained for each kind of field. First of all we must consider resonances of the 
gravity field with the dilaton contribution. We must consider the fluctuation 

G'mn = Gmn + huN, (4.1) 

where Hmn represents the graviton in the axial gauge h^M = 0. We must also assume that 

h^^ = d^h^^' = 0, (4.2) 

because we are interested on the transverse modes. With these considerations we can 
obtain the equation, 

|_e2(^-S) + ^2{A-B)j^,d_^ ^ 2e2(^~^) {X' - A'B' + 2{A'f) - a^j V = 0, (4.3) 

with the definition = r]^'^djj,dy. Performing now the separation of variables /i^jy = 
hf_iu{x)ip{y), we get the equation 

- /(y) + B'ij'iy) + 2 - A'B' + 2{A'f^ ^{y) = m\''^''-^H{y). (4.4) 

It has been shown in [27] that the zero mode is localized. Using now the relation between 
A and B and our previous transformations (3.7) we get the Schroedinger equation 

{-^ + W = mV, (4.5) 

with potential given by 

U=l^^'^{X' + \{^fy (4.6) 

Now we study resonances in a dilaton free background. The strategy here is the same 
as that used to find the background solution without the dilaton contribution. The final 
result is obtained just by setting i? = to arrive at the equation 

|-e2^(^)^ + 2e2^(^) (a" + 2(^')') " 5^} V = 0, (4.7) 
with the same definitions as before. We get therefore the equation 

~ + (2X' + 4(^')') } V'(y) = m2e-2^(^)V(2/). (4.8) 

In order to study the massive spectrum and resonances we must transform the above 
equation into a Schroedinger equation. This is easily obtained with the transformations 
(3.7) and we obtain 

d2 



^ U{z) \ ^{z) = m^ij{z), (4.9) 



with potential U{z) given by 



U{z) = -e^"" ( 2A' + 5(^0' ) • (4.10) 



- 11 - 



1.0 r 



-I.OL 



Figure 6. Potential of the Schrodinger like equation of gravitational field with dilaton (lined) and 
without dilaton (dotted). 



(a) 

log(7-) 


-5 

-10 

-15 

-20 

(b) 

log(r) 


-5 
-10 
-15 
-20 



Figure 7. Logarithm of the transmission coefhcient (a) without dilaton and (b)with dilaton. 

We show in the Fig. 6 the graphics of the gravity potential. As said in the second section, 
this potential has the general behaviour lim2„j.-|-oo U {z) = 0. We must compute transmission 
coefficients to study the possibility of massive gravity with the method of transfer matrix. 
We show in Fig. 7 the result graphically. 

From the Fig. 7(b) we can see a peak of probability for massive gravitons near m? = 
0.02125. Therefore this massive mode is to be expected to interact with the membrane. 
We should point that we have the absence of resonances without dilaton, as we can see 
in Fig. 7(a). This result agree with [17] although the author has used a slightly different 
topological defect. It is important to note that the appearance of resonances is model 
dependent. 
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5 The Scalar Field Resonances 

In this section we study the massive modes of scalar fields. This can be considered the 
simplest application of the method. 

Initially we take the action for the scalar field ^ coupled to gravity, 



^ J d^xdyV^e-^^G^'^ dM^dN^, (5.1) 



2 

where the indexes M, go from to 4. We continue by finding the equations of motion, 

V'^'d^d^^ + e-2^-^+^"aje^^-^-^"9j,$] = 0, (5.2) 

where we have separated the extra dimension from the others. In the next step we use the 
following separation of variables: 

^{x,y) = x{x)i^{y)- (5.3) 

Writing $ as above we arrive at the following equation for the y dependence: 

+ -B'- A7r)-f }V = -m2e2(^-^)^. (5.4) 
dy dy 

Now we must use our relation between B and A and the transformation (3.7) to get our 

Schroedinger equation 

+ f/(z)l^ = mV, (5.5) 



dz^ 

where the potential U (z) assumes the form 

m = e'^'^ ((^ - l^)J^{y? - (f + \)J^'{y)) , 

where a = -15/4 - XVSM^. 

For the dilaton free case we just fix i? = vr = to obtain the equation o motion 

r]'"'df,d^<!> + e-^^^y'>dy[e^^^y^dy^] = o, (5.6) 

and performing again the separation <I>(x,y) = x{x)'^iy) we get for the y dependence: 

{t^ + = -m^e-^^V- (5-7) 

dy dy 

Now using (3.7) we obtain the Schroedinger equation 

^ + U{z)\^ = m^^, (5.8) 



where the potential U (z) assumes the form 



U{z) = e 



2A 



'-liA'f + lA" 



(5.9) 



We show in Figs. 8 and 9 the graphics of the potential and the logarithm of the transmission 
coefficient for the scalar field. As pointed in the last section the presence of the dilaton 
can induce the appearance of resonant peaks. In Fig. 9(c) we see a very rich structure of 
resonances for AVSM^ = 40. This indicates that, for this model, some massive scalar field 
have a peak of probability to interact with the membrane. 
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Figure 8. Potential of the Schrodinger like equation of scalar field scaled by 1/10. (a) with dilaton 
(Hned) for AVSM^ = 2 and without dilaton (dotted), (b) with dilaton for AVSAf^ = 40. 




Figure 9. Logarithm of the transmission coefficient (a) without dilaton, (b) with dilaton for 
AV3Af3 = 2 and (c) with dilaton for XVSIVP = 40. 



6 The Vector Field Resonances 

Here we must consider the next important field, which is the gauge field. It has been show 
that this field is localized only in the Dilatonic background. Anyway, for completeness, we 
give an analyses of the resonances for the dilaton free case. These cases have already been 
studied in [27] [45], but we analyses it here in the light of the Transfer Matrix Method. 

The coupling between the dilaton and the gauge field is inspired in string theory and 
we have for the action 

Sx = I d^'xV^e-^-lYM.M.Y^''^^ (6.1) 

where YmiM2 = d[MiXM2] is the field strength for the 1-form X. The equation of motion 
is: 

dM{V^g''''g'''^e->^^HpQ) = 0. (6.2) 
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We can use gauge freedom to fix Xy = d^X^ = and we are left with the following two 
type of terms 

= ^[mi^M2]' (6-3) 

Yy^ = dyyX,] = dyX, = X'^. (6.4) 
The equation of motion takes the form 

g^^ymM2 + g((-B+A^)5^[g(2A-B-A^)^/M2j ^ (g_5^ 

and performing the same separation of variable we get 

- (-2A' + B' + Att') ^}V^ {y) = -m2e^(^-^)^ (y) . (6.6) 

Now we can get our Schroedinger equation by performing the transformation (3.7) and use 
the relations between vr and A. We arrive at 

U{z) = ((^ - l^)A'{yf - (| + ^)A"(y)) , , 

where a = -7/4 - XV^KP 

Now we repeat the same lines used throughout all the paper to reach the Dilaton Free 
case. First of all the action is given by 

Sx = j d'xV^iYM.M.Y''^^^ (6.7) 

with the same definitions as before. The equation of motion can be obtained from the 
previous case by fixing i? = tt = and we get 

Q^^ymM2 + dy[e^^X'^^] = 0, (6.8) 
and performing the same separation of variable we get 

+ 2A'^}ij (y) = -m'e-^^i^ (y) . (6.9) 



Now we can get the Schroedinger equation by performing the transformation (3.7) to arrive 
at 

U = e^^{^A'^ + U"}. (6.10) 

We show in the Fig. 10 graphics of the vector field potential. In the Fig. 11 we present the 
logarithm of the transmission coefficient. As commented in the introduction, the coupling 
with the dilaton induces a rich structure of resonances. Without the dilaton we see that 
we do not have any resonances. With the inclusion the dilaton we see in in Fig. 11(b) that 
we find one resonance peak very close to zero. When we enlarge the value of the coupling 
constant we get six peaks of resonances, a result which corroborate our claim that very 
massive modes can also interact with the membrane. 
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Figure 10. Potential of the Schroedingcr like equation of vector field scaled by 1/10. (a) With 
dilaton (lined) for XVmP = 2 and without dilaton (dotted), (b) with dilaton for AVSM^ = 40 
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Figure 11. Logarithm of the transmission coefficient, (a) Without dilaton, (b) with dilaton for 
XVUP = 2 and (c) with dilaton for XVmP = 40 

7 The Kalb-Ramond Resonances 

The dilaton coupling here is the same as that for the scalar and vector fields. The reason 
is that the dilaton can be seen as the radius of a compactified extra dimension and all the 
fields get an exponential factor coming from the measure. Therefore we have for the action 

Sx = j d''x^e-''^YM,M,M,Y^'''''-'^\ (7.1) 

where YmiM2M3 = f^[Mi-^M2M3] is the field strength for the 2-form X. The new equation 
of motion is: 

5A./(V=^/^^<7^«/^e-^-i7pQH) = 0. (7.2) 
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We can use gauge freedom to fix X^-^y = d^^X^^^,^ = and we are left with the following 
two type of terms 

^y^^l^^2 = ^[yX^^^^] = dyX^^^^ = X'^^^2- (7.4) 
The equation of motion takes the form 

and we can now separate the y dependence of the field using 

X^'^^'^ {x",y) = B>''^'^ {x'')^{y). (7.6) 
Defining y^i^2 = Y^^^^^ij), where Y stands for the four dimensional field strength, we get 

^^^ymM2 _^ rn'^B^'^^' = (7.7) 
and the differential equation which give us information about the extra dimension, namely 

- iXAy) + B'iy))'^ = -m^e^(-(^)-^(^))^(,). (7.8) 

Using our previous relation between B,tt and A and performing the transformation (3.7) 
we get 

+ [/(z)l^ = mV, (7.9) 



(7.10) 



where the potential U (z) assumes the form, 

9 \ , ^ a 3 



U{z) = 62^ 

where, 



4 64 / ^ ^ \ 2 8 ' 



P=^ + l,a=\-V3M~^X. (7.11) 
Now we introduce the action for the Kalb-Ramond field without the dilaton as 

Sx = j d^xV^lHMNLH^'^''], (7.12) 

and we can again obtain the equations of motion setting B = ir = X = 0. The result is 
given by: 

9^^ymM2M3^(y) _ e2A^^%)^^2M3 = 0. (7.13) 

The function tp{y) carry all information about the extra dimension and obeys the following 
equation: 

^ = -m'e-'^(y^m. (7.14) 
dy^ 
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We must make now the standard transformation (3.7) to get 

|-^+^w}?^(^)=m2v^, (7.15) 

where 

U{z) = e'^[\iA'f + l{A")]. (7.16) 

We show in the Fig. 12 the graphics of the Kalb-Ramond field potential. In Fig. 13 we 
show the logarithm of the transmission coefficient. As in the scalar and vector field, the 
resonance only exist for large coupling constant. Note that, very similar to the vector field 
case, the number of resonances increase very much when consider a large value for the 
coupling constant. 
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Figure 13. Logarithm of the transmission coefficient, (a) Without dilaton, (b) with dilaton for 
AV3A/3 = 2 and (c) with dilaton for An/SM^ = 40. 
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8 The Fermion Resonances 



The next important case to be studied is the fermion field. The issue of zero mode local- 
ization of this field has attracted much attention recently [33-43]. Here we are interested 
only in the study of resonant modes. They depend strongly of the form of the potential 
and, therefore, of the model considered. Here we must study the model similar to that 
found in [27] and [38], therefore we must give only a short review on how to arrive in 
the Schroedinger equation that determines the massive modes. The action considered is 
slightly modified to include the dilaton coupling, as in our previous sections. The action is 
given bellow 

5i/2 = J d^x^e-^- [^T^'Dm^ - r]^F{<P)^] , (8.1) 

where Dm = Qm+^m, being the spin connection. The F*^ = e^^^7^ are the Dirac ma- 
trices in the five dimensional curved space-time and ej^'^ are the vielbeins: e j^^'^ e r]^^ = 
. After some manipulations with the spin connection we get the equations of motion 



[-f^'^^ + {dy + 2dyA) - r/e^F(,/.)] ^' = , (8.2) 

where ^^d^ is the four-dimensional Dirac operator on the brane. It is interesting to note 
here that the dilaton coupling just modifies the equation of motion throughout the B 
factor of the metric. Therefore the value of A will be irrelevant to analyses the resonance 
spectrum. Now, just as in the gravity case we perform the decomposition 

^(X,y) = \Y.i^Ln{x)fLn{y) + Y.'t'Rn{x)fRn{y)\ ■ (8.3) 

\ n n / 

In the above we defined ^^'4'Ln{x) = —tpLn{x) and 'y^tpRn{x) = tpRn{x) as the left and 
right-handed fermion fields respectively. Just as pointed in the second section we obtain 
the equation for a massive field in four dimension 

l^d^tpLnix) = mntpRnix) 

and 

Yd^ijRnix) = mnlpLnix), 

with mass defined by 

[dy + r?e^F(</.)] fLniy) = m^e^'-^fRniy) , (8-4) 

[d, - r?e^F(0)] fRniy) = -m„e^-^/L„(y) . (8.5) 

From the these equations we can obtain the Schroedinger equation for each chirality 

[-dl + UL{z)]fL = m^fL, (8.6) 
[-d! + Ur{z)] Jr = m^fn , (8.7) 
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with potential given by 

Ul{z) = e2^-2^ (r?2e2^F2(0) - rje"" dyF{cl>) - T^e"" {dyA)F{<P)) , (8.8) 
Un{z) = e2^-2^ (r?2e2^F2(0) + r/e^a,F(<A) + r/e^(9^A)F(</.)) , (8.9) 

with dz/dy = e^~^. With this at hand we can now study resonances. As said in all 
previous section, to obtain the case without the dilaton we just fix i? = 0. We show 
in the Fig. 14, Fig. 15, Fig. 16 and in the Fig. 17 the graphics of the fermion potential and 
the logarithm of the transmission coefficient for F{(j)) = (p. We note that if the potential 
looks like the double well barrier, we can observe resonances. In the case of the right 
fermion it looks like a one single barrier and consequently we cannot observe resonances. 
In Fig. 15, with r] = 0.5, we can see that regardless of the dilaton coupling we do not have 
any resonances. In Fig. 16, for = 10 we see that some peaks appear for the left handed 
fermion. Therefore it looks that the resonant structure here depends mostly on the value 
of r]. 
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Figure 14. Potential of the Schroedinger like equation of fermion with 77 = 0.5. (a) Left without 
dilaton (dotted) and with dilaton (lined) and (b) Right without dilaton (dotted) and with dilaton 
(lined). 

9 Conclusions 

In this work we have addressed the issue of resonances in Randall-Sundrum-like models. 
This matter has already been studied before for many authors, but the general prescription 
is to use an non-physical normalization and try to solve the equations of motions. This 
approach could hide the physics involved in the problem and we propose to study it by 
the transfer matrix method of computing transmission coefficients. As pointed in the 
introduction this is a common technique used in condensed matter systems. It is used for 
example to find the peaks of currents in a given semi-conductor. Roughly speaking, the 
method consists in taking the potential of the problem and approximate it by a series of 
potential barriers. Then, for each barrier (for each interval) it is computed interactively 
the transfer matrix, obtaining finally the total transfer matrix. With this method we have 
computed resonances of gravity, fermions and gauge fields. 

For the gravity field, we show in the Fig. 7 the graphic of the logarithm of the trans- 
mission coefficient. From that picture we can see that we have no peak of probability 
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Figure 15. Logarithm of the transmission coefhcient of fermion with 77 — 0.5. (a) Left without 
dilaton, (b) Left with dilaton, (c) Right without dilaton and (d) Right with dilaton. 




Figure 16. Potential of the Schrodinger hke equation of fermion with rj = 10. (a) Left without 
dilaton (dotted) and with dilaton (lined) and (b) Right without dilaton (dotted) and with dilaton 
(lined). 



for massive gravitons in Fig. 7(a). This correspond to dilaton free case. Therefore a mas- 
sive mode is not expected to interact with the membrane. For the dilaton case, Fig. 7(b) 
we find a resonance peak. It is important to note that the appearance of resonances is 
model dependent. For the simplest case, the scalar field, we show in Fig. 9 the graphic 
of the logarithm of the transmission coefficient. The presence of the dilaton can induce 
the appearance of resonant peaks. In Fig. 9(c) we see a very rich structure of resonances 
for A V 3M^ = 40. This indicates that, for this model, some massive scalar field have a 
peak of probability to interact with the membrane. In the case of the vector gauge field, 
the coupling with the dilaton induces a rich structure of resonances. Without the dilaton 
we see that we do not have any resonances. With the inclusion the dilaton we see in in 
Fig. 11(b) that we find one resonance peak very close to zero. When we enlarge the value 
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Figure 17. Logarithm of the transmission coefficient of fcrmion with 77 ~ 10. (a) Left without 
dilaton, (b) Left with dilaton, (c) Right without dilaton and (d) Right with dilaton. 

of the coupling constant we get six peaks of resonances, a result which corroborate our 
claim that very massive modes can also interact with the membrane. The Kalb-Ramond 
field behaves as the scalar and vector fields, i.e., the resonance only exist for large coupling 
constant. We note that, very similar to the vector field case, the number of resonances 
increase very much when consider a large value for the coupling constant. For the case 
of fermionic fields we note that if the potential looks like the double well barrier, we can 
observe resonances. In the case of the right fermion, it looks like a one single barrier and 
consequently we cannot observe resonances. In Fig. 15, with ry = 0.5, we can see that, 
regardless of the dilaton coupling, we do not have any resonances. In Fig. 16, for r] = W 
we see that some peaks appear for the left handed fermion. Therefore it looks that the 
resonant structure here depends mostly on value of rj. 

Furthermore, there are some applications of the models studied here related to AdS/CFT 
correspondence: there are evidences supporting the idea of domain wall/QFT, i.e., a cor- 
respondence between gauged supergravities and quantum field theories in domain walls. 
Important applications are related to the study of quark-gluon plasma via gravity-duals. 
Therefore, we hope to adapt this numerical method to attack such kind of problems. 
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